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SOME EXPERIMENTS AT HIGH SUPERSONIC SPEEDS ON THE 
AERODYNAMIC AND BOUNDARY -LAYER TRANSITION 
CHARACTERISTICS OF HIGH-DRAG BODIES 
OF REVOLUTION 

By Alvin Seiff , Simon C. Sommer, and 
Thomas N. Canning 


SUMMARY 


Measurements have heen made at Mach numbers of It-.O and 8.3 of the 
drag, static stability, damping in pitch, and boundary-layer transition 
characteristics of a number of high-drag bodies of revolution of a type 
that might be used for high-speed entry into the earth’s atmosphere. 

It was found that the static stability of the test models was inherently 
very great with centers of pressure as far aft as 91 percent of the 
length from the nose. The damping coefficients measured for a round- 
nosed 60° cone were found to be adequate for rapid damping of pitching 
oscillations of a full-scale vehicle. Boundary- layer transition on 
these models was found to occur at low Reynolds numbers compared to what 
has been obtained on slender models at the same free-stream conditions. 
Raising the Mach number had a favorable effect on the extent of laminar 
flow, but did not significantly increase the Reynolds number of transi- 
tion when based on local flow properties . The reasons underlying the 
early transition were studied and are discussed. 


INTRODUCTION 


In reference 1, it was shown that bodies of moderate size and weight 
can enter the earth’ s atmosphere at hypersonic velocity with less aero- 
dynamic heat input if the wave drag is large than if the body is slender 
and aerodynamically efficient . This is due to the fact that wave drag 
converts missile kinetic energy to heat in the atmosphere, whereas the 


2 


NACA RM A 56 IO 5 


large skin-friction drag associated with slender "bodies converts missile 
kinetic energy to heat in the missile. This finding focuses interest 
on the hypersonic aerodynamic characteristics of bodies with large wave 
drag. Many things are poorly understood in this flow regime, and the 
authors of this report chose to investigate the following: 

1. Boundary -layer transition characteristics of high-drag bodies. 

2. Static and dynamic flight stability of high-drag bodies. 

3 . Adequacy of Newtonian theory for predicting force and moment 
characteristics of these bodies. 

The importance of boundary-layer transition at hypersonic speeds 
is generally known. Turbulent boundary layers transfer significantly 
larger amounts of heat than laminar boundary layers, the exact compari- 
son depending on particular conditions. Therefore, the difficult problem 
of coping with aerodynamic heating during atmospheric entry is alleviated 
if the boundary layer remains laminar. 

There has been considerable speculation on the likelihood of laminar 
flow on high-drag bodies at hypersonic speeds . This speculation has 
been based on subsonic experience with spheres, supersonic and hyper- 
sonic wind-tunnel experience with slender bodies, supersonic wind-tunnel 
experience with spheres, and laminar stability theory. The subsonic 
experience showed that a strong negative pressure gradient such as 
occurs on a sphere is very favorable for keeping the boundary layer 
laminar. The supersonic experience with slender bodies showed that 
heat transfer into the body significantly increases the stability of 
laminar flow. Spheres tested supersonically at nearly zero heat transfer 
(refs. 2 , 3 , and h) have given indications of laminar boundary layers 
at moderate Reynolds numbers. Laminar stability theory has been used 
(see, e.g., ref. 5 ) to conclude that high-drag pointed cones of large 
included angle will, for all reasonable values of wall temperature in 
flight, have stable laminar boundary layers at hypersonic speeds. Since 
pointed cones may require rounding to avoid loss of the tip due to over- 
heating in flight, the effect of rounding off the cone is of interest. 
This has been discussed in references 6 and 7 and was found to be favor- 
able to laminar flow in one important sense - the local Reynolds numbers 
are appreciably reduced due to rounding. Thus, several indications lead 
to the expectation that laminar flow should not be difficult to retain 
on spheres and round-nosed cones at hypersonic speeds. 

Two factors have appeared to work against keeping the boundary- 
layer laminar on a missile entering the atmosphere. The first is angle 
of attack. The missile can be expected to enter the atmosphere with 
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some angular misalignment, and ■thereafter to undergo a pitching oscil- 
lation. Slender bodies at angle of attack have been shown (ref. 8) to 
experience transition on their sheltered sides at angles of attack as 
small as 1 °. The expectation for low-fineness-ratio bodies, refer- 
ence 8, is that they will not be as sensitive to angle of attack. How- 
ever, the changes in the boundary layer due to angle of attack are so 
complex that this cannot be accepted with assurance. A second factor 
opposing laminar flow is surface roughness. Hypersonic missiles may 
be subjected to dust collisions, surface oxidation, and surface ablation 
in flight, thus producing, in some degree, a roughened surface. It is 
therefore particularly important to know the effect of surface roughness 
on transition for this application. 

With this background information in mind, it was decided to make 
tests of hemispheres, pointed cones with 60° included angle, and 60° 
cones with the nose rounded off spherically. As the tests proceeded, 
additional test configurations were added to investigate particular 
points. The test conditions were Mach numbers of 4 and 8 with the 
models at room temperature and with stagnation air temperatures of 
2200° and 4300° Rankine, thus simulating the flight condition of large 
heat-transfer rate into the body surface. 

Little is known about the static and dynamic flight stability of 
high-drag bodies. For cones of large included angle, the center of 
pressure is known to lie somewhat behind 2/3 the length from the nose 
(the position for slender cones) . For the other test bodies only the 
Newtonian theory was available for estimating the static stability , 
but the accuracy of such estimates was not known. Moreover, no measure- 
ments of dynamic stability derivatives for high-drag bodies at hypersonic 
speeds are known to the authors. Therefore exploratory measurements. of 
the aerodynamic characteristics of the test bodies were made to see if 
they were satisfactory for practical application. 

Work has been done on the drag and pressure distributions of high- 
drag configurations at hypersonic speeds. Reference 9 , for example, 
gives pressure distributions on round-nosed cones and hemispheres, while 
references 2 , 3 , and 4 give similar data for spheres. Drag data has 
been reported for cones of large included angle (ref. 10) and spheres 
(refs. 11 and 12 ). Drag measurements were likewise made in the present 
test and are compared in the report with the preceding data and with 
Newtonian theory. 

Preliminary results of that portion of the present tests dealing 
with transition were previously reported in reference 13. 
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NOTATION 


frontal area, ft 2 
drag coefficient, 

<lco A 


lift coefficient. 


lift 

<lc<A 


dC T 

lift-curve slope, — -, per radian 

da 


pitching-moment coefficient, P atching moment 

< 3 00AI 


dC- 


pitching-moment-curve slope, — per radian 

da 


damping -moment coefficient , 


Sc 


m 




damping-moment coefficient. 


Sc 


m 




normal-force coefficient. 


normal force 

q«A 


normal-force-curve slope, per radian 

da 


pressure coefficient. 


P - Poo 
( loo 


pitot-pressure coefficient. 


Pt - Po 


diameter, ft 


nose diameter, ft 
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e 

f 

h 



Jo 

ki 

2 

k 2 

k 3 

k 4 

Z 

m 

M 

nji 

T 


p 

Pt 


Rn 


R h 


Rs 


Naperian base 
frequency, cps 
roughness height, in. 

transverse moment of inertia about y axis through center-of- 
gravity, mu 2 , slug ft 2 

Bessel function of the first kind 

damping constant defined in equation (5)> ^ 1 


frequency- determining constant defined in equation (AlO) , ft 2 

damping constant defined in equation (B12) 

frequency- determining constant defined in equation (B13) 

reference length, ft 

mass, slugs 

Mach number 

number of observations, np, in which a given station on the body 
is turbulent, compared to the total number of observations, N 

static pressure, lb/ft 2 

total pressure behind normal shock wave, lb/ft 
free- stream dynamic pressure, lb/ft 2 

free-stream Reynolds number based on nose diameter, 

Pi u ih 

Reynolds number based on roughness height, — — — 

Reynolds number based on diameter, 


Poo u oo d n 


p-,Ui s 

Reynolds number based on distance along model surface, — 

P’1 
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Reynolds number based on axial distance from leading edge, 

Pl*U- 2.X 

i^i 

distance along model surface, ft 

distance along model surface to point of transition, ft 
time, sec 
temperature, °R 
velocity, ft/sec 

velocity along flight path, ft/sec 

velocity at which missile enters earth 1 s atmosphere, ft/sec 
coordinate along the model axis or along the flight path, ft 
distance from nose to center of gravity, ft 
distance from nose to center of pressure, ft 
altitude, ft 

angle of attack, radians except where noted 
envelope of pitching oscillation, deg 
constant, l/22,000, ft -1 
flight -path angle, radians 
boundary- layer thickness, ft 

initial boundary -layer thickness at stagnation point, ft 

wave length of pitching oscillation, ft/cycle 

coefficient of viscosity, lb sec/ft 2 

air density, slugs/ft 3 

datum air density, 0.0034 slugs/ft 3 

radius of gyration, ft 

angle between body axis and horizontal, radians 
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0c 

©E 

i 


cone half -angle, deg 

flight-path angle at which missile enters earth’s atmosphere, 
deg 

inclination of body surface to free-stream direction, deg 



Subscripts 

00 free-stream conditions 

1 conditions at outer edge of boundary layer 

w conditions at model surface 

i initial conditions 

R resultant 


Derivative Notation 

(*),(") first and second derivatives with respect to time, respectively 

(’),(") first and second derivatives with respect to distance, 
respectively 


APPARATUS AND TEST TECHNIQUE 
Facilities 


The test facilities used for this investigation were the Ames super- 
sonic free-f light wind tunnel (ref. l 4 ) and the Ames supersonic free- 
f light underground range. In the former, gun-launched models are flown 
upstream through the test section of a supersonic wind tunnel to produce 
hypersonic test conditions. The wind-tunnel test section is equipped 
with four shadowgraph stations at intervals along the flight path 

to record the motion history of the models. Reynolds number variation is 
obtained by variation of the supply air pressure. The underground range 
is a more conventional ballistic-type firing range in which models are 
fired through still air at atmospheric pressure. The instrumentation 
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consists of seven shadowgraph stations at 11-foot intervals. In both 
facilities, test models were launched from a 1.75~i n ch bore diameter 
smooth-bore gun. 


Models 


The basic model, the 60° included-angle cone with a spherical tip, 
had a tip radius equal to one -third the cone-base radius, and a base radius 
of O.863 inch. The geometry behind the base varied, depending upon the 
launching technique used. One type had a short, threaded base extension, 
with diameter equal to 71 percent of the base diameter, which screwed 
into a split nylon sabot. After leaving the gun, the sabot halves sepa- 
rated leaving the model in free flight. The second type had a slender 
shaft extending from the center of the base (diameter about 8-l/2 percent 
of the base diameter and O.75 inch long) which not only supported the 
model in the sabot but also aided in the determination of the angle of 
attack from the shadowgraphs. Both types were made of 7075 - T6 aluminum, 
and are shown seated on their respective sabots in figure l(a) . A third 
variation of this model was geometrically identical to the second type; 
however, the forward 90 percent of the model length was constructed of 
magnesium and the last 10 percent plus the shaft was constructed of 
high-strength steel in order to move the center-of-gravity position aft. 

The auxiliary models used are represented in figures l(b) and l(c) . 
Figure l(b) shows a 60° included-angle pointed cone and a hemisphere 
with the same base diameter as the round-nosed cone of figure l(a) , and 
also shows a 3/8-inch-diameter sphere. The cone and hemisphere were 
made of 7O75-T6 aluminum, and the sphere was a commercial steel ball 
bearing. A slight modification of the hemisphere, a spherical segment 
subtending a total angle of 214°, was also employed. Illustrated in 
figure l(c) is a series of aluminum models that were tested primarily 
to investigate boundary-layer transition. The model shown at the top 
of the figure is a segment of a pointed ogive. The next model differs 
only in that it has a rounded nose, the rounded portion (enlarged and 
illustrated to the left of the model in the figure) being a section of 
a sine wave of revolution faired into the basic ogive with continuous 
first and second derivatives at the tangent point. The model shown at 
the bottom of the figure is a spherical-tipped 30° included-angle cone. 
This model, unlike the others, was launched using the vacuum technique 
described in reference 14 for holding the model on the sabot. 


Surface Finish 


The methods of surface preparation and inspection used in this 
investigation were patterned after and developed from the methods 
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described in reference 8. Three classes of surface finish were employed. 
These will he described, and then the techniques used to inspect them 
will be outlined. 

Type I .- Successively finer grades of emery polishing paper were 
used, with care exercised at each step to remove the polishing scratches 
of the previous step. For the final polish the finest available grade, 
designated 4/0, was used. This procedure resulted in a finish with cir- 
cumferential scratches as much as 10 to 20 microinches deep (as deter- 
mined with a surface interferometer to be described below) . The scratch 
width or spacing was determined from photomicrographs to be about ^>00 
microinches or 0.0005 inch. Scratches of this dimension are perceptible 
to the naked eye, as can be seen in figure l(a) , which shows models with 
Type I surface finish. 

Type II.- Models with Type I finish were repolished with maximum care 
using grades 3/0 and k/0 paper in an effort to improve the uniformity of 
the surface. Scratches were comparable to those in the Type I surface. 

Type III .- Models with Type II surface were further polished with 
aluminum- oxide abrasive of 20-microinch particle size. The circumferen- 
tial scratches were reduced to a depth of 3 to 5 microinches and were 
normally less than 0.01 inch long. However, the Type III surfaces had 
more three-dimensional roughness, usually pits left by pulling out alloy- 
ing inclusions, than Types I and II. Figure l(d) shows a model with a 
Type III surface. 

One of the techniques used to inspect the model surfaces was direct 
examination with a metallurgical microscope. A photomicrograph at 500X i 
typical of either a Type I or a Type II finish is shown in figure 2(a). 

The polishing scratches are very conspicuous. A Type III surface at the 
same magnification is shown in figure 2(b) . The scratches are now greatly 
reduced but still perceptible. The alloying inclusions and pits referred 
to above show up as dark spots in the photomicrograph. 

To measure the height of the roughness, a surface interferometer was 
employed. It consists of a partly silvered optical flat placed in con- 
tact with the model surface, with the region near the point of contact 
in view through a microscope. Monochromatic light is used for illumi- 
nating the mirror from above. Fringes are then observed in the microscope. 
Photomicrographs of the fringes are shown in figures 2(c), 2(d), and 2(e). 
The fringes result from interference between light reflected from the 
mirror and light reflected from the model surface. They can be shown to 
be contour lines of constant distance between the model surface and the 
mirror. Adjacent fringes mark contours which differ in elevation by half 

■^he bread, dark line running across the middle of the picture is a 
wire shadow used for measurement of rough surfaces. It has no significance 
in the present case. 
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the wave length of the light. In the present case, green light from the 
mercury-arc spectrum was used, and the half wave length is approximately 
10 microinches. Where the contour lines are smooth, the surface is 
smooth. Where the contour line is displaced, for example, by half a 
fringe space, there exists a 5-microinch irregularity in the surface. 

Figure 2(c) shows an interferogram of a Type I surface. The worst 
scratches and ridges are seen to be about 10 microinches high and, in 
a few places, the depth may exceed 10 microinches. The Type III surface, 
shown in figure 2(d) , has roughness nowhere greater than about 5 micro- 
inches. The root -mean- square roughness would, of course, be considerably 
less. The stagnation point of a model with Type III surface is shown in 
figure 2(e). A scratch ranging from 3 to 10 microinches deep runs directly 
across the stagnation point. Examination of representative pictures like 
these three led to the selection of peak roughness amplitudes given 
earlier . 


Range of Test Variables 


When the models were fired through still air in either the wind w 

tunnel or underground range, the Mach number was limited to about 4.5. 

Tests in the wind tunnel with an air- stream Mach number of 2 gave test 
Mach numbers up to 9 .3* The model velocities ranged from 3500 to 6000 
feet per second and were controlled by varying the propellent powder 
charge. The nominal test Mach numbers, temperature conditions, and 
Reynolds numbers based on free- stream properties and maximum diameter 
are listed below for each of the models tested. 


Moo 

Model type 

Roo*10~ 6 

Free-stream 
static air 
temperature. 
Too, °R 

Stagnation 
temperature 
(computed) , 
°R 

Tw/Too 

3 

60° cone 

3.0 

530 

1450 

1.0 

4 

Round-nosed 60° cone 

Hemisphere 

Pointed ogive 

Sine -wave -rounded ogive 

Round-nosed 30° cone 

3.8 

530 

2200 

1.0 

8.3 

Round-nosed 60° cone 
60° cone 
Hemisphere 

4.8 

295 

4300 

1.8 

9.3 

3/8-in. -diameter sphere 

1.1 

295 

5300 

1.8 
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A pertinent variable for transition tests, the ratio of local wall 
temperature to local static air temperature just outside the boundary- 
layer, T w /T x , is omitted from the above table because T x varies along 
the body, except in the case of the pointed cone, so that no single value 
of T w /T x will describe conditions for the entire body. However, it can 
be noted that T x is always higher than for the models tested, 

and therefore that T w /T x is always less than the values tabulated above 
for V-V 


Data Reduction 


The metnods used for measuring force and moment characteristics will 
he described in this section. Methods used to detect the location of 
boundary-layer transition will, because of their intimate relation to the 
results obtained, be included in the later section where boundary- layer 
transition is treated in detail. 

Aerodynamic forces and moments studied included drag, static sta- 
bility, dynamic stability, and (briefly) lift-curve slope. These charac- 
teristics were determined from the shadowgraph records of model motion 
as a function of time and distance. Drag coefficients were computed from 
the deceleration of the models by the procedure given in reference l4. 
Static stability was computed from the frequency of pitching oscillation. 
To determine the frequency, a damped sine wave. 


a 


- 

= e 2 ( C x cos 2it ^ - C 2 sin 2jt ^ 


( 1 ) 


is fitted by least squares to the available observations of a as a 
function of distance traveled, x. An example set of data and the fitted 
sine wave are given in figure 3 • From the pitching frequency, the static 
stability derivative, Cn^, was computed by use of the relation 

(2rtf) 2 I y k x 2 Iy + CmqClciftxA 1 + GpC^ftpAIy 

q^AZ 2p oQ AZ 2m 2Zm 2 



(2a) 


which is derived in Appendix A. This equation assumes linear pitching 
moment and damping proportional to pitching rate, but it differs from 
equations that have been given previously in that the effect of model 
deceleration on the pitching motion is accounted for. Including this 
effect leads to the result that the pitching frequency is proportional 
to velocity, but that the wave length of oscillation. A, is constant. 
There is, then, a significant advantage in using distance as the 
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independent variable in equation (l) rather than time. Furthermore a 
small-order term involving Cp is found to occur in equation (2a). This 
term, however, and the two preceding terms are found to be negligible for 
the present test conditions (less than l/2 percent of the first term, 
total) even though the velocity loss in flight was relatively large. 

The resulting equation used was 


Ci% - 


_ q^Al J 


8 it 2 ly 

Poo Al 


(2b) 


where, in the case of the countercurrent air stream, A must be defined 
as the wave length relative to the free-stream air. 

The pitching-moment data obtained in this way were broken down to 
give normal-force-curve slopes and center-of -pressure locations by use 
of the familiar relation 


Cm a ~ C N a 


x cp ~ x cg 
l 


(3) 


Since both center-of -pressure location and normal- force -curve slope are 
unknown, the data from two models with different center -of -gravity loca- 
tions are required. (As was noted in the section on models, a configu- 
ration using heavy metal in the base section was employed to give a rear 
c.g. location.) In one case, however, it was not necessary to use two 
models to obtain this breakdown of the pitching-moment data. By chance, 
a model flew with sufficient pitching amplitude, 12°, to cause the center 
of gravity to fly a measurably swerving course. From the curvature of 
the flight path, the lift-curve slope, Cx^, was computed by the method 
of reference 15. The value of Cjj was then obtained from the relation 

c N a = ^ + C D (4) 


which can be shown to hold accurately for small angles of attack from 
consideration of the vector diagram of aerodynamic forces. The measure- 
ment of Cn^ for this round by the methods outlined in the preceding 
paragraph, and the application of equation (3) for center of pressure 
gave the complete static-stability data from one round fired. 

The measurement of dynamic stability goes back to the determination 
of the constant, kj_, in equation (l) by obtaining the best fit of a 
damped sine wave to the angle-of-attack data. In Appendix A, it is shown 
that ki is related to the dimensionless aerodynamic damping coefficients 
as follows: 
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Hence can be readily computed once k 1? Cp, and are 

known. Of even more interest, however, is the fact that the group of 
coefficients in the bracket of equation (5), designated £ for brevity, 
is the same group obtained by Allen in reference l 6 . According to this 
reference, £ must be negative to give a well-damped motion of a missile 
on re-entry. It must likewise be negative to give a convergent oscilla- 
tion in the present test. Hence (except for possible differences in 
radius of gyration between test and flight), convergence of the motion 
in the test indicates convergence of the motion for re-entry. The rate 
of convergence, however, is not in general the same, since it is governed 
by multiplying factors involving air density, model frontal area and 
mass, etc., which are not the same either in algebraic form or in magni- 
tude for range tests and re-entry flight. Hence, the combined group of 
dimensionless aerodynamic coefficients, £ , represents the most general 
way of stating the dynamic characteristics of a configuration. 


RESULTS AND DISCUSSION 


The results to be presented will concern the force and moment char- 
acteristics and the boundary-layer transition characteristics of the test 
configurations. The former data can be presented and discussed briefly 
and in a straightforward manner because there are few points of interpre- 
tation involved. Consequently this will be done first. Following this, 
the boundary-layer transition data will be presented. As will be seen, 
this section will include a description of the shadowgraph evidence used 
to detect transition, a discussion of various factors possibly influencing 
the results, and some speculation as to what the test results might 
indicate for objects of larger scale. 


Force and Moment Data 


The force and moment measurements were made not so much with the 
idea of cataloging the characteristics of a wide range of configurations 
as with the idea of providing a few spot checks of the theoretical methods 
available. Thus, drag measurements were limited to four configurations 
at two Mach numbers; normal force and center of pressure, to two 

2 The similarity of form results from including the effect of decel- 
eration in both cases. (Note the appearance in eq. (5) of Cp.) The 
effect of deceleration is to reduce the dynamic pressure, and therefore 
it is destabilizing with respect to pitching amplitude. 
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configurations at two Mach numbers; and damping- in-pitch measurements 
were limited to the basic configuration at one Mach number. These 
measurements, however, may indicate what can be expected from bodies 
of this category, and also the ability of the available theory to 
predict the individual characteristics. 

Drag.- Total-drag coefficients of the round-nosed 60° cone, pointed 
60° cone, hemisphere, and sphere were measured at Mach numbers near k 
and 8 and are plotted as a function of Mach number in figure k. The 
round-nosed 60° cone and the pointed 60° cone have drag coefficients that 
are essentially equal. The drag coefficients of the hemispheres and 
spheres are substantially higher and nearly independent of Mach number 
over the test Mach number range. 

A comparison of the experimental total-drag coefficients with theory 
is shown in figure 5(a) for the round-nosed cone, hemisphere, and sphere. 
Two methods of calculating total-drag coefficients were used; the method 
of reference 17, and a modification of the Newtonian impact theory of 
reference 18. The method of reference 17 consists of adding the wave 
drag of the spherical tip, from experiment, to the theoretical wave drag 
of the cone (from ref. 19) and subtracting the wave drag of the conical 
tip. Additional allowance is made for base drag and skin-friction drag. 
The modification of Newtonian impact theory consists of substituting the 
pitot-pressure coefficient, for the coefficient 2 in the expression 

for pressure coefficient as a function of flow deflection angle, 0 S . 

Thus we use 


Cp ~ Cp^sin 2 0 s (6) 

instead of 

Cp = 2 sin 2 0 s (6a) 


This modification tends to give more nearly the correct pressure coeffi- 
cients on the blunt portions of the body and becomes exact at the stag- 
nation point of round-nosed bodies. Again allowance is made for base 
drag and skin-friction drag. Precise estimation of base drag and skin- 
friction drag is not required, however, since the major contribution to 
the total drag of these bodies is wave drag. The base pressure was there- 
fore arbitrarily taken to be 0.3 of the free- stream static pressure, and 
a skin-friction coefficient of 0.0025 based on wetted area was arbitrarily 
selected. From these assumptions, the base drag coefficient at M^ =8.0 
becomes 0.0l6, and the skin-friction drag coefficient based on frontal 
area becomes 0.005, compared to an estimated wave-drag coefficient of 
O.562 for the round-nosed cone and 0.915 f° r the sphere. 
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The calculations by the method of reference 17 give the drag of the 
round-nosed 60° cone within 5 percent. The drag estimated using modified 
Newtonian theory is too low for the round-nosed 60° cones by as much as 
10 percent , and too high for the spheres by about 5 percent, although the 
agreement in the case of the former model seems to improve as Mach number 
is increased. The use of unmodified Newtonian theory, equation 6(a), 
gives estimates about 10 percent higher, and thus improves agreement with 
experiment for the round-nosed cone at the lower Mach number, but worsens 
agreement in the case of the sphere. 

The total-drag data for the sharp 60° cones of this test are plotted 
in figure 5(t>) along with the experimental drag curve for the 60 ° cone- 
cylinders of reference 10, which agree with the present measurements very 
satisfactorily. Of course, the exact conical-flow theory of reference 19 
provides precise prediction of the wave drag for this model (ref. 10) , 
but to test further the modified Newtonian theory, equation (6) was 
applied to this case. The result (including friction and base drag) is 
shown in the figure by the solid curve. It underestimates the measured 
drag by as much as 18 percent. Thus, in three comparisons, the modified 
Newtonian theory is found both to underestimate and overestimate drag, 
with a maximum discrepancy of 18 percent. Unmodified Newtonian theory 
is subject to errors of similar magnitude. The method of reference 17 
appears somewhat more accurate for the limited cases to which it can be 
applied. 

Static stability .- No presentation is made here of the values of 
Cn^ measured, since they would apply only for the particular center-of- 
gravity locations used. Instead, the quantities, C ^ and x cp , which 

specify the static stability for any desired center-of-gravity location 
are given. 

Experimental values of Cjj for the round-nosed 60° cone are shown 

in figure 6 and compared with theory. The normal-force-curve slopes fall 
somewhat below the value indicated by slender-body theory, 2.0, but this 
still gives the order of magnitude, even under these conditions of extreme 
departure from slenderness. The modified Newtonian theory tends to fall 
somewhat below the data as shown. Unmodified Newtonian theory would agree 
somewhat better. The theory for a pointed 60° cone (ref. 19) falls in the 
middle of the data for the two Mach numbers. To what extent this could be 
expected to continue for cones of different included angle and increasing 
amounts of bluntness is unknown. 

Center-of -pressure data for the round-nosed 60° cone and the pointed 
60 ° cone are given as a function of Mach number in figure 7* (Note that 
in the case of the pointed cone, square symbol, the theoretical value of 
CN a vas use( 3- with the measurement of to obtain the result shown.) 

In both cases, the center of pressure falls approximately 90 percent of 
the length from the nose. Prediction of the center of pressure of the 
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round-nosed cone by Newtonian theory is surprisingly accurate, within 
3 percent of the length, in view of the discrepancies in prediction of 
drag. The result for the pointed cone is in very satisfactory alignment 
with conical flow theory, and provides confirmation for the following 
center-of -pressure equation for pointed cones with attached shock waves: 

= | (l + tan 2 0 c ) 

In the case of the hemispheres, no static stability measurements 
were made. However, it is of interest to mention that the center-of- 
pressure forces on the front face will fall at the center of curvature - 
that is, at the base. The bodies were designed with this expectation in 
mind, and their flight characteristics were very satisfactory, a proper 
nose-forward attitude being consistently maintained. 

The above results indicate that there should be little difficulty in 
obtaining statically stable flight with bodies of this type at high speeds. 

Dynamic stability .- Damping- in-pitch information was obtained from 
two round- nosed 60° cones at a Mach number of 1. Round 1 was fired in 
the wind tunnel (air off) and gave four angle -of -attack observations in 
one cycle of oscillation. Round 2, which is the one shown in figure 3> 
was fired in the underground range, and provided seven observations of 
angle-of-attack spread over 3*5 cycles. The data from these two rounds 
are summarized in the following table: 


Round 

no. 

Mach 

no. 

x cg / z > 
percent 

Amplitude 
ratio per 
cycle 

c mq + 

(exp.) 

c m q 

(theory) 

(exp.) 

(theory) 

1 

3.98 

70 

0.78 

- 0.70 

- 0.23 

-5.1 

- 2.0 

2 

3.91 

70 

• 7^ 

-.85 

-.23 

-6.2 

- 2.0 


Considering the differences in the distribution and extent of the raw 
data, and the minimal definition of the oscillations, the agreement of 
the two results is fairly good. The agreement of Newtonian theory with 
experiment is not very satisfactory. The theory gives an estimate of 
Cjj, only, which is about one third as great as the measured values of 


Cm q + Cm d * 


The predicted value of £ is likewise about one third as 

great as measured. The accurate prediction of £ will be very difficult 
because £ involves sums and differences of the aerodynamic coefficients 
and precision depends on extreme accuracy (or compensating errors) in the 
individual coefficients. Further investigation will be required to 
determine fully the adequacy of the theory. The measured coefficient, 

C m + C^, was negative, indicating a tendency to damp the motion. The 
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combined coefficient , £ , was negative as well because the tendency for 
deceleration to destabilize the motion (Cp in eq. ( 5 )) was, in this 
case, fully canceled by the damping due to plunging (C^ in eq. (5))* 

At the time these measurements were made, it was not clear that the 
coefficients obtained would or would not lead to a well-damped re-entry 
flight. Therefore, the calculations described in Appendix B were under- 
taken to consider the flight of an example missile. Friedrich and Dore 
(ref. 20) had obtained equations for the pitching motion during re-entry 
of a body with zero aerodynamic damping and zero damping due to plunging. 
They considered only the effect of the variation in dynamic pressure, 
which damps the motion at the beginning but causes divergence to occur 
at low altitudes when q begins to decrease because of velocity loss. 

The envelope of a motion of this type is shown in figure 8 by the curve 
labeled "dynamic-pressure effect - aerodynamic damping = 0." Approximate 
equations developed in Appendix B were then applied to determine the con- 
tribution of aerodynamic damping, using experimental values of the damping 
coefficients. This led to the amplitude variation indicated in the figure 
as "contribution of aerodynamic damping." The resultant motion is then 
obtained by multiplying the ordinates of the two curves which leads to the 
curve labeled "total effect." The aerodynamic damping coefficients 
measured are found to be very beneficial at altitudes below 100,000 feet. 
They are, in fact, sufficient to overcome the tendency for divergence at 
low altitude. 

Subsequent to these calculations, the linear equations of planar 
pitching motion during re-entry were solved with only minor approximations 
as reported in reference l6. Calculation of the motion envelope by the 
methods of that reference led to identically the same "total effect" 
curve as that given in figure 8. (See Appendix B.) 


Boundary-Layer Transition 


Transition to turbulence of the boundary layers on the test models 
was studied from the spark shadowgraphs. Because the boundary layers and 
disturbed-flow fields were thin, and because light refraction effects 
were strong, transition locations were not so immediately apparent in 
these shadowgraphs as in those of the slender bodies of reference 8. It 
was possible in most cases, however, to find evidence of the type of 
boundary- layer flow. 

Evidence of turbulence.- It was found that boundary- layer turbulence 
on these bodies could produce one or more of at least four indications in 
the shadowgraph: 

1. In some cases, the turbulence could be directly observed as in 
reference 8. This required a relatively thick boundary layer in a region 
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free of light distortion from the bow shock wave and was therefore usually 
limited to positions toward the back of the test bodies. 

2. Where the turbulent boundary layer was very thin, the turbulent 
eddies frequently caused minute light filaments to be refracted into the 
model shadow, producing the impression of hairiness at the model edge. 

3. Sometimes Mach waves produced by the turbulence could be seen 
well outside the boundary layer when the turbulence itself was pot evident. 
The turbulence waves are of greater slope than body-fixed Mach waves. 
Sometimes they occur as isolated single waves in the case of bursts of 
turbulence. Other times, when transition is relatively steady, a field 

of turbulence waves will occur and it will be bounded by an envelope cor- 
responding to a body-fixed Mach wave originating at the transition point. 

4. The condition of the boundary layer behind the model base can 
frequently be determined from its appearance. The laminar layer gives a 
thin, smooth, steady line which can persist far behind the base, as in 
figures 9(a) and 9(d). The turbulent layer gives an indefinite, notice- 
ably turbulent line which diffuses rapidly, as in figures 9(l) and 9(c). 

These features are all indicated schematically in the accompanying 
sketch. In the course of the discussion, shadowgraph pictures from the 
test showing some of these characteristics will be introduced. Turbulence 
was considered to be present in a given case if any one of the above 
indications could be observed. 
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Round-nosed 60° cone ,- Shadowgraph pictures typical of those obtained 
with the round-nosed 60° cones with Type I surface finish at a Mach number 
of 8.3 and a Reynolds number of 4.8 million, based on maximum diameter 
and free- stream air properties, are shown in figure 10. These pictures 
contain examples of two of the turbulence indications cited above. In 
figure 10(a) the shadowgraph sensitivity is low and, as a result, turbu- 
lence in the boundary layer is directly visible in the original negative 
as, for example, on the top surface starting at about midlength. In 
figure 10(b) the view of the surface is obscured because of higher sensi- 
tivity but a laminar line can be seen behind the base on the bottom 
meridian. Thus, it can be concluded that there is fully laminar flow on 
the bottom but not on the top. 

This asymmetry in transition was noted for two models tested at this 
condition. One side was laminar throughout the flight, independent of 
angle of attack and regardless of whether the laminar side was at the 
moment windward or sheltered. On the other side, the boundary- layer flow 
off the base was consistently turbulent. This was taken as evidence of 
an influence of roughness on transition, since roughness was known to be 
almost certainly asymmetric (see, e.g., fig. 2(e)), whereas no other 
significant characteristic of the model or air stream was known to be 
persistently asymmetric. 

Graphical presentation of the data from the model of figure 10 and 
one other model at a Mach number of 8.3 is given in figure 11. In 
part (a), the location of transition is plotted against angle of attack. 
Three kinds of symbols are employed. The circles without arrows indicate 
direct observation of transition to turbulence on the model or behind the 
base. The circles with downward- facing arrows indicate that transition 
occurred at some unknown location ahead of the base. The x f s indicate 
the observation of turbulence waves, usually single waves from transient 
bursts of turbulence. The figure indicates that transition occurred on 
the model on a majority of the meridians observed. The effect of angle 
of attack appears to have been relatively moderate compared to that 
observed with the slender models of reference 8. Transition apparently 
began to move forward on both windward and sheltered sides at angles of 
attack above 2°. 

The local Reynolds numbers at transition are presented in fig- 
ure ll(b) . For transition on the body, the Reynolds number was based 
on surface length to the transition point and air properties just outside 
the boundary layer at the transition point. To compute the air properties 
the total pressure outside the boundary layer was assumed to be that 
behind a normal shock wave. The static-pressure distribution on the 
spherical portion was computed by modified Newtonian theory, equation (6a) 
The static-pressure distribution on the conical portion was approximated 
by smooth fairing between the terminal pressure of the sphere and the 
static pressure a.t the base, assumed to be that of a pointed cone. Even 
fairly large changes in the assumed pressure distribution would not 


20 


NACA RM A56I05 


significantly alter the Reynolds numbers shown in figure ll(b) . Behind 
the base, pressure was assumed equal to l/lO free-stream static pressure, 
and this led to a value of local Reynolds number per unit length behind 
the base about l/lOO of that in the free stream. Thus, the sizable runs 
of laminar flow behind the base shown in figure ll(a) represent only 
slight increases in transition Reynolds number over the value for tran- 
sition at the base. The local transition Reynolds numbers are of the 
order of 1 million or less. 

Additional tests of this model configuration were made at a lower 
Mach number, near 4.0, and at a slightly lower Reynolds number, about 
3.8 million based on diameter and free-stream conditions. Surface 
finishes of Type I (four models) and Type III (one model) were employed. 

A shadowgraph print of this model at Moo = 4.0 is shown in figure 12 . 3 
Notice that the background of tunnel-wall turbulence is absent for this 
case, since the test was made in still air, and there is a consequent 
gain in shadowgraph clarity. Also the lower Mach number increases the 
separation of the shock wave from the surface (which is somewhat contrary 
to the frequently expressed idea that Mach number has no influence on 
these flows). As a result, turbulence can be directly observed wherever 
it occurs, except on the spherical nose. Two other things were noted in 
the original negative. First, the hairy condition was seen where the 
boundary layer is thin, and this develops into visible turbulence as the 
boundary- layer thickness increases. Second, the turbulence behind the 
model base could be clearly seen both by looking at the edge of the base 
flow and by looking through the base flow. 

Transition data collected from such pictures are presented in fig- 
ure 13 . Part (a) is comparable to figure ll(a) except that the angle- 
of -attack range extends up to 11° instead of 3°« It also differs in that 
transition was found to occur forward on the body in every case observed. 
The effect of angle of attack again appears to be small. Separating the 
data according to meridians gave a somewhat more systematic effect of 
angle of attack. The two or three observed points on any individual 
meridian showed in almost every case that increasing the angle of attack 
tended to move transition forward. Differences in roughness of the indi- 
vidual meridians caused them to have differences in transition location 
which are largely responsible for the vertical spread of the data in 
figure 13 (a). 

The data for the model with the Type III finish are shown by the 
filled symbols. No improvement in the extent of laminar flow was pro- 
duced by the additional polishing. This may be due to the roughness 
remaining in the Type III finish after even the most careful polishing 


3 The bright lens-shaped image located about l/2 inch behind and 
parallel to the model base is due to light from the shadowgraph spark 
reflecting off the model base. 
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(see figs. 2(d) and 2 (c)). The roughness of the Type III finish tends 
to be three-dimensional , and hence may be as destabilizing to the laminar 
boundary layer as the scratches of the Type I finish at this test 
condition. 

The local transition Reynolds numbers represented by the data of 
figure 13(a) were computed by the method previously described and are 
shown in figure 13(b) • Tbe most interesting point to be noted is that 
while the extent of laminar flow on this configuration was less at a 
Mach number of 4 than at a Mach number of 8, the local Reynolds number 
of transition was nearly the same at both Mach numbers, about 1 million. 

Pointed 60° cone .- For comparison with the above data, limited tests 
were made of a pointed 60° cone with Type I surface finish at a Mach 
number of 8.3. The Reynolds number based on free-stream properties and 
diameter was 4.8 million, the same for the pointed cone as for the 
round-nosed cone at this Mach number. 

A shadowgraph from this model test is shown in figure 14 (a) . Close 
examination of the model edge did not reveal any of the indications of 
turbulence cited earlier. In particular, the edge of the model appeared 
very clean and smooth, free of visible eddies or hairiness. Calculations 
were made to estimate the thickness of a turbulent boundary layer on this 
cone at the test conditions. The thickness computed to occur at the model 
base was 0.020 inch for transition at the tip and 0.017 inch for transi- 
tion at midlength, thicknesses which should be easily visible in the 
pictures. It was therefore concluded that the absence of any visible 
turbulent boundary layer and the absence of any burst waves indicated 
that the boundary layer on the surface was fully laminar. Study of the 
wake boundaries did not, however, reveal any laminar line as distinct 
and persistent as in figure 10 (b) . A faint laminar line appeared inter- 
mittently, both near the base and farther back, which suggested that 
transition to turbulence occurred unsteadily along the boundary of the 
separated region. 

To determine the appearance of these flow details in a case of defi- 
nitely turbulent boundary layer, a cone was roughened by sandblasting 
with number 60 carborundum grit. One of the shadowgraphs of this cone 
in flight is reproduced in figure 14 (b) . Three differences from the pic- 
tures of the smooth cone were noted: Turbulence was visible at the model 

edge. No fine line whatsoever was visible at the edge of the base region. 
The minimum section of the wake was advanced from a mean position of about 
1.4 diameters behind the base for the polished model to about 0.8 diameter 
behind the base for the rough model (see also fig. 9). These differences 
tended to support the conclusion that the smooth cone had fully lami nar 
boundary layer. 

The local slant-length Reynolds number of the cone at this flight 
condition was computed to be 6.7 million. The local Reynolds number of 
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transition would therefore he something in excess, of 6.7 million. The 
angle-of-attack range of the data was from 0.9° to 3*8°, hut the data 
were too limited to permit the evaluation of possible angle-of-attack 
effects . 

Comparison of results from the rounded cone and the pointed cone 
showed that, for the particular conditions of these tests, rounding the 
cone reduced the average extent of laminar flow. This occurred even 
though the local Reynolds numbers were appreciably higher in the case of 
the pointed cone. The surface finishes of the two configurations were 
nominally the same. Nevertheless, the round-nosed cone gave indication 
of roughness effect in that transition differed on different meridians, 
whereas the pointed cone with Type I finish showed no tendency for rough- 
ness to cause transition. It is therefore indicated that the rounding 
increased the sensitivity of the laminar boundary layer to surface 
roughness as a contributing cause of transition. 

Hemisphere .- Hemispheres were used as a second comparison configu- 
ration, thus giving the case of no rounding (pointed 60° cone), inter- 
mediate rounding (round-nosed 60° cone), and full rounding (hemisphere) 
to be compared. Tests of hemispheres with a Type I surface finish were 
made at the two basic Mach numbers, 4 and 8.3, and at free- stream Reynolds 
numbers based on diameter comparable to those of the round-nosed cones. 

The shadowgraph pictures of the hemispheres at the higher Mach number 
were very difficult to interpret. One of them is shown in figure 15(a). 
The boundary layer on the surface was not clearly visible, but in the 
heavily exposed pictures, the model edge was seen to be hairy, to within 
7° of the stagnation point in one picture, corresponding to a transition 
Reynolds number of about 20,000. The appearance of the boundary layer 
behind the model base was somewhat puzzling. It seemed to have laminar 
characteristics for a short distance (0.1 or 0.2 of an inch) before 
becoming distinctly turbulent. This may indicate an effect of strong 
expansion (along the sphere) of the type observed and discussed by 
Sternberg in reference 21 for the expansion at the corner of a 60° cone- 
cylinder. The relative turbulence will be decreased by rapid expansion 
if the turbulent velocity fluctuations grow, as would be expected, at a 
smaller rate than the velocities in and at the edge of the boundary layer. 
If this is the correct explanation for the apparent laminar flow just 
behind the model base, it suggests the possibility that laminar or inter- 
mediate heat-transfer characteristics might be occurring toward the rear 
part of the hemisphere. This suggests an interesting topic for further 
research, a test of the heat transfer on a sphere with turbulence induced 
near the stagnation point. 

For contrast with the above model test, a hemisphere roughened by 
sandblasting was tested, and shadowgraphs like figure 15(b) were obtained. 
While the hairy condition at the model edge was somewhat intensified, and 
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the length of laminar line behind the base was somewhat reduced, the main 
impression to be derived from comparing these pictures with those of the 
smoother model was that they were substantially the same. 

Tests at the lower Mach number of b of the 214° spherical segment 
model (fig. 15 (c)) gave definite evidence (hairiness and visible turbu- 
lence) of turbulence on the front face at locations as far forward as 
within 40° of the stagnation point. The local Reynolds number at this 
point was 0.8 million. Because of optical distortion, it could not be 
determined how far ahead of this point transition took place. Toward 
the back of the model, the apparent intensity of turbulence decreased, 
but the flow coming off the model base remained in this case definitely 
turbulent . 

Comparison with earlier data from slender-model tests .- At this 
point, it was apparent that transiton was occurring at far lower Reynolds 
numbers on the round-nosed models than on slender pointed models tested 
under similar free-stream conditions (refs. 8 and 22). The slender models 
had been able to retain laminar boundary layers with appreciable two- 
dimensional roughness present - orders of magnitude greater than that on 
the present models. The reasons for the relatively poor performance of 
the present round-nosed models in retaining laminar boundary layer were 
naturally the subject of considerable speculation. As was discussed in 
reference 13 , one aspect of the flow on nonslender bodies of this type 
which seemed significant was the compression to surface Mach numbers much 
lower than free-stream Mach number. From the viewpoint of laminar sta- 
bility theory (ref. 23 ), the large regions of subsonic flow on the round- 
nosed bodies could not be completely stabilized by cooling, even if the 
disturbances were restricted to the more favorable, two-dimensional kind. 
Experimentally, it had been shown (ref. 22) that at low wall temperature, 
reducing the Mach number at the boundary- layer edge increased the sensi- 
tivity of the laminar boundary layer to roughness- induced transition. 

For these reasons, attention was focused on the low surface Mach numbers 
occurring on the test bodies as a primary factor causing early transition. 
Another thing that aroused some interest was the thickness of the boundary 
layers on these bodies. The occurrence of roughness effects on relatively 
smooth surfaces seemed consistent with the expectation that the boundary 
layers would be very thin. Therefore, calculations of boundary- layer 
thickness were undertaken. Finally, unusual aspects of the flow over 
these bodies were considered to see if some additional destabilizing 
influences could be found. Several such influences will be discussed 
subsequently. 

Because it was suspected that the low values of local Mach number 
were instrumental in causing early transition, it was desired to test a 
body with higher local Mach numbers, but still with a fineness ratio 
comparable to that of the other bodies tested. This led to the tests 
described below. 
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Pointed ogives ,- A model was designed to have a minimum surface 
Mach number of 2 at a free-stream Mach number of 4. This, of course, 
required that it be pointed. Continuous expansion along the surface to 
a Mach number of 2.3 at the base was provided. The resulting model was 
the pointed ogive shown in figure l(c). Two of these were prepared, one 
with Type I finish and the other with Type II. Since the supersonic 
free-f light wind tunnel was by then shut down for modifications, the 
tests of this model were made in the underground range and were limited 
to a free-stream Mach number of 4. The Reynolds number based on diameter 
was the same as in the earlier tests at this Mach number. 

The results obtained are shown in figures 1 6 and 17* Figure 16(a) 
shows the model with Type I surface and shows bursts of turbulence (as 
indicated by the burst waves) occurring very far forward. Figure 16(b) 
shows the model with the Type II surface. It is free of turbulent burst 
waves, but shows very fine-grained turbulence beginning at about one fifth 
the length from the nose. (This latter picture was obtained at reduced 
shadowgraph sensitivity to reduce the distortion due to the bow shock 
wave. The method used to reduce the sensitivity was to bring the film 
very close to the plane of flight, within about 2 inches.) The data from 
the two rounds are presented graphically in figure 17« Since no angle- 
of-attack effects were noted, the distribution of observed transition 
points is presented according to the method of references 8 and 24. The 
number of cases in which a given body station is found turbulent, Up, 
divided by the total number of observations made, N, is plotted as a 
function of body station. The mean location of transition is seen to 
have been relatively near to the tip where the relative accuracy of deter- 
mining the transition location is poor. The difference shown in the 
mean location of transition on the two models is therefore not believed 
to be significant. 

The reasons why these models failed to retain laminar flow were not 
clear at the time but, subsequently, a correlation of the roughness neces- 
sary to cause transition as a function of local Mach number showed that 
this test was still in the region where roughness controls transition. 

This point will be further discussed in the section on "Comparison of 
Roughness Height to Boundary -Layer Thickness." 

Round-nosed ogive .- To obtain additional information on the effects 
on transition due to rounding a pointed body, a round nose was applied 
to the pointed ogive. The nose diameter was selected to give a low 
Reynolds number layer thicker than the boundary layer, according to the 
argument of reference 6 . A sine-wave rounding was employed to provide 
continuous first and second derivatives at the juncture. A shadowgraph 
of the round-nosed ogive is shown in figure 18 , and the transition loca- 
tions read are compared in figure 19 with those of the pointed ogive. 

The effect of rounding, in this instance, was to increase the extent of 
laminar flow as predicted in reference 6 . It is to be noted, however, 
that in the earlier case of the 60° cone at a Mach number of 8 , rounding 
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reduced the extent of laminar flow. Thus, depending on particular fac- 
tors, rounding may or may not increase the extent of laminar flow. It 
appears that effects of rounding, such as changes in local Mach number 
and pressure distribution, can be as important as the effect on local 
Reynolds number and must therefore be considered. 

Comparison of the pointed and round-nosed ogives on the basis of 
local transition Reynolds number is made in figure 20 , which is like 
figure 19 except that the body stations are now represented by their 
local Reynolds numbers. This presentation inverts the relation between 
the pointed and round-nosed ogives - the round-nosed models show slightly 
lower average Reynolds number of transition. This indicates the presence 
of some moderately adverse effects of rounding. 

A comparison of round-nosed ogives with Type II and Type III surfaces 
is given in the figure. The two surfaces gave almost the same transition 
pattern. Rather than indicating no effect of roughness, this indicates 
(more generally) that the two surfaces are of equivalent roughness. 

Round-nosed 30 ° cone .- In order to investigate a case intermediate 
in fineness ratio between the preceding bodies and the very slender bodies 
of references 8 and 22, the spherical-nosed 30° included- angle cone of 
figure 1 (c) was tested in the underground range. The Mach number was 4.3 > 
the Reynolds number based on free- stream properties and diameter was 
4 million, and the surface finish was Type II. Shadowgraphs are presented 
in figure 21. 

Transition to turbulence still occurred relatively early. This can 
be seen in the shadowgraphs both from visible turbulence and burst waves. 
The data read from the shadowgraphs are presented in figure 22 , where 
observations from the four individual meridians observed are kept sepa- 
rate and denoted right, left, top, and bottom. Apparently, some system- 
atic differences between the transition on individual meridians occurred, 
although these differences were not so pronounced as in the case of the 
round-nosed 60 ° cone. In the present case, for example, the performance 
at zero angle of attack was best on the right-hand meridian, worse by a 
factor of 3.5 on the left, and intermediate on top and bottom, possibly 
because of microscopic differences in smoothness of these meridians. 

The effect of angle of attack shown in this figure seems to be 
generally similar to that of figures 11 and 13, namely, a moderate effect, 
with transition tending to move forward on both windward and sheltered 
sides as angle is increased. It is thought that in the case of more 
extensive laminar flow, a stronger effect of angle of attack might be 
observed on this body. 

The local Mach number occurring on this body just ahead of the base 
was estimated to be near 2.0 from the slopes of Mach waves visible in 
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the shadowgraph. The local Reynolds numbers at transition, based on 
wetted length and estimated air properties outside the boundary layer, 
were again of the order of 1 million. 

Comparison of Roughness Height to Boundary-Layer Thickness 

The roughness-height data described in the section on surface rough- 
ness are not of themselves sufficient to describe the roughness conditions 
of these tests, but must be considered in relation to the boundary- layer 
thicknesses on the test models. In reference 13, a roughness parameter 
to relate roughness height, h, to laminar boundary- layer thickness, 5, 
for cases of distributed surface roughness of uniform height on a flat 
plate was proposed. Since this parameter was only briefly mentioned in 
reference 13, it will be further treated in this section, with the objec- 
tive of applying it to the correlation and interpretation of the present 
test results. 

Pointed bodies .- In the case of a pointed cone or sharp-edged flat 
plate with uniform surface roughness, a difficulty arises in attempting 
to specify a characteristic value of h/5, since this ratio will vary 
from infinity at the leading edge to values generally less than 1 near 
the base. Therefore, the specification of h/5 must either be arbitrary, 
as in references 8 and 25, or in the form of a graph giving the distribu- 
tion of h/5. For flat plates and cones, however, it is possible to 
identify a parameter which specifies the relationship of roughness to 
laminar boundary-layer thickness for the entire surface. To find this 
parameter, first note that on the flat plate. 


where C is constant over the surface of the plate , but is a function of 
the boundary- layer-edge Mach number and the ratio of wall temperature to 
boundary- layer-edge temperature. In order to make it possible to use the 
theoretical results of Van Driest (ref. 26) for the constant C, the 
boundary-layer thickness is defined throughout this report as the distance 
from the surface at which u/ui = 0.995* From equation (7)> then, it 
follows that 


so that h/5 is a function of the Reynolds number per unit length, 
Pi u i/Pii an< 3- distance from the leading edge, x. Equation (8) can be 



(7) 



(8) 


written 
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h _ r pi^ih/ii! 1 — - 

5 Lc(p 1 u 1 x/p 1 )J' / x 

(9) 

or 

5 V x C 

(10) 


The right side of equation (10) is seen to be constant for flat plates 
of uniform roughness and a given Reynolds number per unit length. By 
the same token, (h/&) -jR-x is seen to be a function of the Reynolds num- 
ber per unit length and the roughness amplitude, but is constant over 
the surface. 

For cones, equation (10) must be modified by the Hantzche and Wendt 
factor, "\/"3, to account for the reduced rate of boundary- layer growth on 
cones. It becomes 

= ( 11 ) 

where C is evaluated from flat-plate calculations for the Mach number 
and static air temperature at the boundary- layer edge. The Reynolds 
number per unit length is also evaluated at the boundary- layer edge. 

For curved bodies, (h/S) ^R^ must, for precision, be computed from 
estimates of 5 which take into account the effect of pressure gradient. 
The roughness parameter will now depart somewhat from a constant value. 

The variation, however, will not be large for moderately curved, pointed 
three-dimensional bodies of any fineness ratio, which do not depart too 
far from conical shape. The ogive of the present investigation was such 
a case. To compute (h/&) sJRq for this body, the boundary- layer thickness 
distribution was first calculated by the method of Cohen and Reshotko 
(see refs. 27 and 28 and Appendix C) . The Reynolds numbers, R s , based 
on local properties outside the boundary layer and surface length were 
calculated using shock-expansion theory. The roughness height, h, was 
taken to be 10 microinches. Combining these data gave the roughness- 
parameter distribution shown in figure 23 . As anticipated, a nearly 
constant value of (h/b) sT^s is obtained. 

Correlation of critical-roughness data for pointed bodies .- Use will 
now be made of the above roughness parameter for studying the roughness 
to boundary-layer -thickness aspects of the present and other tests. 
Critical values of the roughness parameter will be defined as the values 
above which roughness begins to control the Reynolds number of transition. 
For present purposes, if roughness caused transition to occur at Reynolds 
numbers of the order of 1 to 2 million, it was considered that the rough- 
ness parameter was higher than critical. If laminar flow persisted to 
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Reynolds numbers of the order of 10 to 30 million, it was considered that 
the roughness parameter was below the critical value. Based on these 
definitions, data from the present test and other sources were collected 
to define the variation of the critical roughness parameter with Mach 
number shown in figure 2b for wall-temperature ratios, T v /T ± , of the 
order of 1. The three points at Mach numbers of 3*5* ^.9, and 6.9 were 
obtained by cross-plotting data from figure 8 of reference 13. The points 
selected as critical were the points at which an increase in roughness 
started to move transition forward. The two upright triangles are data 
from the present investigation. The pointed 60° cone at a free- stream 
Mach number of 8.3 bad a surface Mach number of 3*0 and fully laminar 
flow to a local length Reynolds number of about 7*10 6 , so the surface is 
considered smoother than critical (arrow up). The pointed ogive at a 
free- stream Mach number of k had a local Mach number of 2 ahead of tran- 
sition and showed early transition, so the surface is considered rougher 
than critical (arrow down) . The cluster of data at Mach numbers between 
1 and 3 represents a broad range of test conditions but, nevertheless, 
demonstrates considerable consistency with the above data. The data of 
Sternberg (ref. 29) and that of Low (ref. 5 ) each give a point on the 
critical roughness curve which is in good alignment with the data for 
higher Mach numbers. Additional points from the data of Sternberg (at 
a Mach number of 2.7), and Rumsey and Lee, reference ^0 , (at a Mach num- 
ber of 2.9) represent conditions of extensive laminar flow and fall in 
the region of smoother-than-critical surfaces. The data of Czarnecki 
and Sinclair, reference 31, and Czarnecki, Robinson, and Hilton, refer- 
ence 32, are again for smoother than critical surfaces at local Mach 
numbers near 1.5, and indicate a somewhat higher position for the criti- 
cal roughness curve than would be deduced from the other data. However, 
the large differences in the conditions of these various tests with 
respect to roughness types, pressure gradient, wall temperature ratios, 
and facilities do not lead to the expectation of quantitative agreement. 
Finally, the data obtained by Luther (ref. 33) in investigating the 
amount of roughness needed to trip the boundary layer near a model tip 
show, as would be expected, that the roughness required to bring about 
immediate transition at a Reynolds number of 9b , 000 is somewhat greater 
than that required to influence transition initially. Nevertheless, the 
trend of variation with Mach number in his data is very similar to that 
shown by the solid curve. It appears that the critical roughness parame- 
ter increases rapidly with increasing local Mach number at supersonic Mach 
numbers above 2. 

The significance of figure 2b in relation to the present test results 
is evident. Although the free-stream test Mach numbers were high, the 
local surface Mach numbers were always less than 3 and- usually less 
than 2. This put the tests into a region where surface roughness effects 
were difficult to avoid. Figure 2b indicates the desirability of obtain- 
ing a more precise definition of the critical values of the roughness 
parameter in the low surface Mach number region where low-fineness-ratio 
missiles are required to operate. 
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Round-nosed bodies .- A calculation of the boundary- layer thickness 
development on the bodies with hemispherical noses (by the method of 
refs. 27 and 28 as discussed in Appendix C) showed that within 30° angu- 
lar displacement from the stagnation point, the boundary- layer thicknesses 
were not represented even approximately by equation (7). In fact, the 
boundary layer on a blunt body has finite thickness initially, and from 
the stagnation point to the sonic point 6 is very nearly constant, as 
is indicated by the calculated thicknesses in figure 25 . This suggested 
that the parameter h/& S p would satisfactorily represent roughness 
effects in this region, and this parameter was adopted for use with the 
round-nosed bodies. 

The calculated ratios of nominal roughness height (taken to be 
10 microinches 4 ) to boundary- layer thickness on the stagnation point are 
tabulated below for the hemispheres and round- nosed cones: 


Body 

Moo 

d n> 

in. 

Poo u oo d n/ Poo> 
million 

5 sp, 

in. 

h/8 S p 

Hemisphere 

3.8 

8.0 

1.73 

1.65 

3.9 

4.5 

0.00121 

.00187 

0.008 

.005 

Round-nosed 60° 

4.0 

.57 

1.3 

.00072 

.014 

cone 

8.3 

.57 

1.6 

.00103 

.010 

Round-nosed 30° 
cone 

4.3 

.25 

.6 

.00048 

.010 


The values of h/b S p computed are of the order of 1 percent, and, judging 
from the observed asymmetry of the boundary- layer transition, roughness 
of this amplitude was sufficient to influence transition. 

Effect of scale .- The maintenance of equal h/b S p provides a basis 
(which is not as yet supported by experiments) for defining round-nosed 
bodies of equivalent roughness but different physical size and Reynolds 
number. The computation of an equivalent roughness does not, however, 
give any information as to the effect on transition Reynolds number of 
increasing the flight Reynolds number at a constant roughness condition. 
This information can be provided thus far only by tests or flights at 
the desired Reynolds number. Nevertheless, it is interesting and informa- 
tive to compute equivalent surface roughnesses to those of the present 
tests for cases of larger scale bodies and increased Reynolds numbers. 

4 The stagnation points of these models were the hardest places to 
polish smooth. All the models examined had several three-dimensional 
pits of the order of 10 microinches deep, within a radius of 0.01 inch 
of the stagnation point. 
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Laminar boundary- layer thickness is, in general, directly propor- 
tional to physical size, and inversely proportional to the square root 
of Reynolds number. Any characteristic dimension may be used as a measure 
of the size, and any single Reynolds number may be used to identify the 
distribution of local Reynolds numbers that exist on a body of particular 
shape. In the following discussion, the nose diameter, d n , will be used 
to specify size, and the free- stream Reynolds number based on nose diame- 
ter, R n , will be used along with the body shape to specify the distribution 
of Reynolds numbers. Then, for given, fixed conditions of Mach number 
and wall-temperature ratio, 

^-<f) 

and 

h _ f _hA N/^n 

5 sp VW f(s/z) sp 

For two bodies of the same shape (but different size and Reynolds number) 
to have the same value of h/5 S p then 


( 12 ) 

( 12 ) 


h •T^n 


"h s/Rn 

dn f ( s/l ) sp 

1 

dn f ( s/l ) cp 


(14) 


where the subscripts 1 and 2 identify the two conditions of scale. If 
Mach number and wall-temperature ratio are the same in cases 1 and 2, 
then 



(15) 


Thus, to make a numerical comparison, consider a body 20 inches in 
nose diameter (60 inches in base diameter) having the shape of the round- 
nosed 60° cone, in flight at M^ = 8 with a free- stream Reynolds number 
based on nose diameter of 15 million, and find the roughness of this body 
which is equivalent to the 10 microinches used in the present test. 

A -J 15x10® = 1 ^/^ 6 \/ 1.6x10 s 
20 0.57 


h = 115 microinches » 0.0001 inch 
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That is, 0.0001-inch roughness on the large-scale body provides equal 
h/§sp to that experienced in the present test. 

In figure 25 , the effect of increasing Mach number on initial 
boundary- layer thickness on round-nosed bodies is shown. The calculated 
initial thickness (ignoring dissociation) for M = 20 is about 3.5 times 
that for M = 3 . 8 . Hence, one might expect that the smoothness require- 
ment on round-nosed bodies might be relaxed by increasing Mach number as 
well as by increasing scale. For the conditions of scale used in the 
above example, the roughness equivalent at M = 20 to that of the present 
test at M = 8 would become 0.00024 inch. 

The situation for scaling roughness on sharp-edged flat plates and 
pointed cones differs fundamentally from that for round-nosed bodies in 
that the flat plate has no characteristic dimension other than the rough- 
ness height (provided that transition occurs within the length of the 
plate) whereas the round-nosed body has as a characteristic dimension, 
the diameter of the nose, on which the boundary -layer thickness depends. 
Thus, for the round-nosed bodies, the scaling relationship suggested 
(eq. ( 15 )) depends on roughness height, nose diameter, and Reynolds num- 
ber; while the roughness parameter expression for flat plates (eq. (lO)) 
depends only on roughness height and Reynolds number per -unit length. 

It would be expected that for flat plates, cones, and perhaps approximately 
for other pointed bodies, that (h/&) n /R q should be kept constant to give 
equivalent roughness, and that therefore 


is the roughness scaling law. 

When it is desired to compare the roughnesses of cones of different 
included angles, or to correlate cone and flat-plate data, then it is 
essential that the Reynolds numbers be based on air properties at the 
edge of the boundary layer. If it is desired only to compare cones of 
identical geometry but of different size and Reynolds number, then free- 
stream properties may be used, since the ratio of free-stream to local air 
properties will be the same if the Mach number is the same. Thus, to 
convert the roughness of the pointed 60 ° cone from the scale tested to 
the same frontal area and flight condition as the preceding example 
(60-inch base diameter and 45 million Reynolds number based on free-stream 
air properties and base diameter) 


R h x - R h 2 



h = 37 microinches 
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Thus, the relaxation of the smoothness requirement with increase in scale 
may be less rapid for pointed bodies than for round-nosed bodies. 

These scaling methods have been discussed at some length because the 
question of proper ways of scaling roughness is a fundamental one for 
studying the transition of boundary layers on large-scale objects from 
small-scale tests. It is emphasized that the procedures used above are 
purely theoretical and have not been investigated experimentally. 


Other Causes of Boundary-Layer Instability 


Attention has been given above to surface roughness as a contributing 
cause of transition in these tests, and it appears that roughness had a 
considerable influence on the results obtained. Other factors, however, 
undoubtedly contributed to causing transition, and these will now be 
discussed. 

Pressure gradients.- Pressure distribution is well known to have 
important effects on boundary-layer transition. This subject has received 
considerable attention in the literature (see refs. 8, 22, 25, 3^> and- 
refs. noted therein). It is therefore necessary to consider the pressure 
fields about the present models. 

The round-nosed 60° cones of the present tests are similar to the 
80° included- angle round-nosed cone of reference 9. In figure 26, taken 
from reference 9, it may be seen that there is a minimum pressure just 
aft of the spherical nose, followed by a fairly quick recovery toward the 
pressure predicted for a sharp 80° cone. Unpublished data from the Ames 
low-density wind tunnel and the data of figure 26 indicate that as the 
test Reynolds number increases, the pressure minimum and subsequent 
recovery become more pronounced. The present tests were made at still 
higher Reynolds numbers so that considerable pressure rise may be presumed. 

Additional evidence of pressure rise along the conical portion of 
the round-nosed 60^ cone was obtained from the shadowgraphs , It was noted 
that the bow- shock -wave profiles were peculiar in that they showed a 
region of reverse curvature - concave outward. This can be seen in fig- 
ures 10, 12, and 21. The reverse curvature was appreciable, and could 
not have been caused by optical distortion since the contour of the bow- 
wave leading edge is recorded by undisturbed light passing through the 
free- stream air just ahead of the shock wave. In the region of reverse 
curvature, the wave slope is increasing, and the wave strength must like- 
wise be increasing. The increase in wave strength is attributed to com- 
pression disturbances feeding up from the body surface, and is therefore 
believed to reflect the condition of rising pressure on the body surface. 
The forwardmost influence of the conical surface on the bow wave has been 
traced in some shadowgraphs of rough models from preliminary tests by 
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following Mach lines attached to "body-fixed roughness. At the intersec- 
tion on the how wave of a Mach wave from the body tangent line, the bow- 
wave profile immediately became straight. Very shortly thereafter, the 
concave profile began. This indicates that the rising pressure on the 
body surface began almost immediately behind the tangent line. 

For purposes of estimating the pressure rise quantitatively, Newtonian 
theory was of no value, since it fails to predict the occurrence of either 
the overexpansion or the pressure rise. If centrifugal force effects are 
included by use of one of the highly approximate procedures available, 
the overexpansion can be introduced, but it will be followed by a dis- 
continuous pressure rise and constant pressure on the conical surface. 

As a more reliable alternative, use was made of sphere pressure distri- 
butions from reference 3 and conical pressure data from reference 19 to 
estimate the pressure rise. Application of this data to the case shown 
in figure 2 6 would, for example, indicate a minimum value of Cp at the 
tangent point of 0.39 Cp max J a va 3 - ue that appears to be in good alignment 

with the pressure data of figure 2 6 for the conical part of the model at 
the highest Reynolds number. When this method of estimation is applied 
to the present round-nosed 60° cone at Moo = 3.8, it is found that the 
pressure coefficient at the tangent line (based on free-stream dynamic 
pressure) is lower by 0.16 than the pressure to which the flow will ulti- 
mately recover on the cone. Converting this result to the form Ap/p min 
(the pressure rise along the surface divided by the minimum pressure along 
the surface) yields a value of 0 . 32 . According to reference 22 (which was 
concerned with slender bodies), this should be enough to assure transi- 
tion if the local Mach numbers are supersonic. It is not known that this 
full pressure rise occurred within the length of model, or to what extent 
the pressure rise that did occur influenced transition. As was noted 
earlier, the hemisphere and ogive models with continually favorable 
pressure gradients had transition Reynolds numbers at M^ = k comparable 
to the round-nosed cones. 

Similar estimates were made of the pressure rise on the round-nosed 
60° cone at M^ = 8 by extrapolating the sphere pressure coefficients of 
Stine and Wanlass to this Mach number. A possible pressure-rise coeffi- 
cient of about 0.17 above minimum was calculated. It is interesting that 
this value is just below the critical value of Ap/p m j_ n given in refer- 
ence 22 and that, for the test at this condition, laminar flow to the 
base was obtained on some meridians. Calculation of the pressures on 
the round-nosed 30° cone at M^ = 4.3 indicated that the terminal pres- 
sure on the spherical nose was slightly greater than that of the pointed 
cone, which would lead to the expectation of favorable pressure gradient 
downstream of the tangent point. This, however, is not consistent with 
the fact that reverse curvature of the bow shock wave was observed, which 
indicates that some pressure rise must have occurred. It further indicates 
that more experimental and theoretical work can be profitably devoted to 
investigating flow fields about bodies of this type. 
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It is suspected that the unfavorable pressure distributions of the 
round-nosed cones contributed to the early transition observed, but there 
is also evidence from comparison with results from the other test con- 
figurations that other causes of laminar instability were present and 
probably of overriding importance. 

Shear -layer effects .- In addition to its effect on pressure distri- 
bution, the flow field outside the boundary layer may also assume impor- 
tance in another way which will now be discussed. Behind the convex bow 
wave lies a region characterized by extreme variations in total pressure 
from one streamline to the next. The region in which these large varia- 
tions occur is quite nicely defined in the shadowgraphs of the round-nosed 
60° cone. As was pointed out above, the juncture of the first Mach wave 
from the cone with the bow wave appears to mark the end of the convex 
portion of the bow wave. The streamline which passes this point lies on 
a sudden change in density gradient normal to the local flow direction. 
Since the shadowgraph is sensitive to changes in density gradient, this 
streamline becomes a visible line in the shadowgraph and shows up clearly 
in figure 10(a). Outside this line the total -pressure variations are 
relatively gradual as evidenced by the comparative straightness of the 
bow wave. The inner region is of interest because of the velocity pro- 
files which can be expected to occur there. Where the total pressure 
varies rapidly across the streamlines, a region of nearly constant static 
pressure is characterized by large velocity gradients and consequently 
by the existence of appreciable shear. This leads to the term "shear 
layer," which will be used in the remainder of this discussion. Layers 
of this kind have been discussed by Moeckel (ref. 6 ) and Ferri and Libby 
(ref. 35 ) . A velocity profile in a shear layer is shown in the follcwing 
sketch as it would appear in inviscid flow, and is qualitatively like 
that between two parallel jets of differing velocities. Wear the nose 
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of such a body in viscous flows , the boundary layer superimposes its 
velocity distribution onto that of the preceding sketch. The resulting 
velocity profile might then be like that sketched below. In this case, 
little interaction between boundary layer and shear layer would be antici- 
pated provided both remained laminar. However, it should be noted that 
this velocity distribution contains two inflection points, one where the 
boundary layer meets the shear layer and one in the shear layer itself. 

It is expected that this velocity distribution would be unfavorable to 
continued laminar flow. 



It is important to note that transition to turbulent flow can occur 
in the shear layer well outside the boundary layer. This phenomenon may 
be seen clearly in the shadowgraphs of figure 21 of the present paper, 
and the occurrence of this transition is indicated in figure 22. Another 
feature of this type of flow is illustrated in figure 27 wherein the 
models are open tubes. The flow on the outside is essentially two- 
dimensional. There is a difference between the two models in the leading- 
edge thickness. In the case of the thin-leading-edge model, figure 27(a*) , 
the boundary- layer growth after transition is relatively slow compared 
to that for the thick-leading- edge model, figure 27 (b). The distance 
from the transition point has been marked off in inches for convenience 
in comparing boundary- layer thicknesses at corresponding stations. In 
the latter case a shear layer (visible in the shadowgraph) has been pro- 
duced by the thick leading edge, and the turbulence quickly envelopes 
this layer, thereby producing a thicker than normal boundary layer. This 
is probably beneficial in that the heat transfer and skin friction of the 
thick turbulent boundary layer should be lower than for a turbulent bound- 
ary layer of normal thickness. Abnormal thickening of the turbulent 
boundary layer was also observed on the original negatives of the round- 
nosed 60° cone. The turbulent boundary layer in these pictures apparently 
envelopes the complete shear layer almost Immediately after transition. 
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Air-stream turbulence .- In the present tests, as in all investiga- 
tions of houndary-layer transition, the effects of air-stream turbulence 
should be considered. In the case of the data for M^ = 4 the tests were 
conducted in still air. The Mach number 8 data on the other hand were 
obtained with a counter current air stream of unknown, but large, turbu- 
lence. When the turbulence velocities are compared to the sum of the 
air and model velocities, the "percent turbulence" is only about one 
fourth as great as when compared to the air velocity alone . However , 
in the subsonic region behind the detached bow wave of round-nosed models 
the compression yields a fivefold increase in turbulence level. This can 
be seen by considering the flow through a normal shock wave where, for 
continuity 


Poo^oo ~ Pl U l 

By the normal- shock relations for 7 = 1.4, 

Pi = 6M^ 

Poo Moo 2 + 5 


Hence, at = 8 



§7 = 5-57 


The resulting increase in percent turbulence probably contributes to the 
difficulty in obtaining laminar flow under these test conditions. 

Gortler instability .- In reference 36 , H. Gortler has treated the 
case of laminar flow over a concave wall and finds this flow configura- 
tion to be favorable to the production of vortices lying along stream- 
lines. In a lecture on laminar boundary- layer instability at Stanford 
University in March 1956, Gortler stated that this phenomenon is not 
limited to boundary- layer flows. Rather, the conditions necessary for 
it to occur are: (a) streamline curvature and (b) velocity gradient 

positive toward the center of curvature. As pointed out by Gortler 
these conditions are met near the stagnation point of a blunt body. 
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In the following sketch, the approximate boundaries of the region of vor- 
tex amplification are shown. (The shock-wave shapes and positions for 
this sketch were obtained from the shadowgraphs of the present test.) 



The shaded area represents the region in question. The outermost stream- 
line involved is probably the one that passes through the point at which 
the maximum flow deflection occurs in the bow wave and the sketch is 
based on this assumption. It can be seen that streamlines inside this 
point will have initially concave profiles. On the body, the curvature 
is always convex, so the region of concave flow may touch the surface 
only at the stagnation point. Therefore it would appear that a thin 
boundary layer would not lie in the region of vortex amplification, and 
that any destabilization would have to occur before the air reached the 
boundary layer or by other indirect means. 

No evidence from the present test confirms or disproves the importance 
of this instability. 


Reynolds Numbers for High-Speed Entry 
Into the Earth 1 s Atmosphere 


If the Reynolds number of a full-scale missile is kept low, more and 
perhaps all of the boundary -layer flow may be kept laminar. It is there- 
fore pertinent not only to calculate the magnitude of the Reynolds numbers 
to be expected in the entry flights of full-scale missiles, but also to 
see what can be done to keep them to a minimum. 

The variation of velocity with altitude of an entering missile is 
given approximately by 
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V = V E e 


_i CpA Pp 
2 m p sin 0 E 


-3y 


for an atmosphere with exponentially varying density (ref. l) . Thus, at 
any altitude, for a particular selection of V E , 0 E , and C^A/m, V is 
known. Furthermore, the density, viscosity, and speed of sound at that 
altitude may he obtained from the assumed atmosphere or from tables of 
air properties as a function of altitude. From these data, the Mach 
number and Reynolds number are readily computed as a function of altitude. 

The velocity and angle of entry, V E and 0 E , are fairly well fixed 
by the flight range of the re-entering missile. However, C-^k/m is a 
parameter under some control by the designer. The effect of this param- 
eter on the Reynolds numbers per unit length is shown as a function of 
Mach number during re-entry in figure 28 for an entrance velocity of 
23,000 feet per second and an entrance angle of 3 O 0 from the horizontal. 
For all the cases shown, peak Reynolds number occurs at the time when 
the Mach number is about 9 . Making the missile light in relation to its 
drag reduces the Reynolds number per unit length. This is accomplished 
either by reducing the weight, increasing the drag coefficient, or 
increasing the frontal area. There is this additional reason, then (in 
addition to those advanced in ref. l) , for using the maximum practical 
value of CjjA/m. This fact was noted by the authors of reference 37 . 


CONCLUSIONS 


From the present tests and analysis of high-drag bodies of revolu- 
tion at Mach numbers near 4 and 8 , Reynolds numbers near 4 million based 
on diameter, and at stagnation temperatures of 2200° R and V 3 OO 0 R, 
respectively, the following conclusions are drawn: 

1. The over-all pressure drag of the bodies tested is predicted 
accurately to within 20 percent in the worst case by modified Newtonian 
impact theory, slightly better than this by unmodified Newtonian theory, 
and more accurately still by the component method of NACA RM A55L21 
where it is applicable. None of these methods, however, is adequate for 
predicting the details of the pressure distribution such as is required 
for interpreting boundary-layer transition data. The component method 
indicates the existence of a region of pressure rise just behind the 
spherical portion of the round-nosed cone, consistent with experiment, 
but does not yield values of the pressure gradient. Impact theory fails 
to predict the existence of the pressure rise. 

2. The static stability of low-fineness-ratio bodies like those 
tested is inherently very great. Center-of-pressure position of the 
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60° cone was at about 0.90 of the cone length from the tip, consistent 
with the theory for low- fineness -ratio cones. Center of pressure of the 
rounded cone was about 0.91 of the length from the nose, as predicted by 
impact theory. 

3. The measured dynamic stability of the round-nosed cone was found, 
from application to an example atmospheric entry of a large-scale body, 

to provide an effective increment to the damping, particularly at low 
altitudes. It was, in fact, sufficient to overcome the tendency for the 
pitching oscillation to diverge as a result of the decrease in dynamic 
pressure at low altitudes. 

4. Boundary- layer transition on the round-nosed bodies was found to 
occur at low Reynolds numbers compared to those that have been obtained 
on slender bodies . The Reynolds numbers of transition based on local 
flow properties were about 1 million at both free- stream Mach numbers in 
spite of the fact that the roughness amplitude was from 0.5 to 2 percent 
of the calculated boundary-layer thickness. Transition at this value of 
the local Reynolds number gave about twice the extent of laminar flow at 
M^, = 8 as at = 4. 

5. For many of the models, transition was observed to be consist- 
ently asymmetric, differing on the different meridians of the model. 

Since surface roughness was also asymmetric, it is considered that sur- 
face roughness influenced transition. Consideration of possible scaling 
relationships indicated that the test models were representative of 
relatively smooth full-scale missiles. 

6. From the available data on the effects of roughness on boundary- 
layer transition on pointed bodies, critical values of a roughness 
parameter were obtained. It was found that as the local Mach number 
increased above 2 the permissible values of the roughness parameter 
increased rapidly. 

7 . Other potentially important causes of instability which may have 
contributed to early transition of the boundary layers are: 

(a) Extensive regions of subsonic flow which are not susceptible 
to complete stabilization, regardless of wall temperature. 

(b) Regions of rising pressure along the direction of surface 
streamlines on some of the test bodies. 

(c) Amplification of the free-stream turbulence by deceleration 
of the main stream velocity through the bow shock wave. 

(d) Instability of the Gortler type in the vicinity of the 
stagnation point. 
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(e) Instability of the shear layer produced behind the highly curved 
bow shock wave. 

8. Transition of the shear layer was noted in several pictures. 

In addition, when transition began in the boundary layer, the shear 
layer was observed to be quickly enveloped to form a thicker-than-normal 
turbulent boundary layer. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Sept. 5 > 1956 
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APPENDIX A 

EQUATIONS OF PITCHING MOTION FOR A MISSILE IN 
DECELERATING FREE FLIGHT AT CONSTANT ALTITUDE 


When a missile decelerates, in flight at constant altitude or in a 
testing range, it experiences a loss in dynamic pressure which tends to 
amplify the oscillating motion and reduce the pitching frequency. There- 
fore, the usual equations for flight at constant speed do not accurately 
describe the motion nor can they, in general, be used to accurately reduce 
observed motions to aerodynamic coefficients. To derive equations which 
include the effect of deceleration, we start with the following differ- 
ential equation of motion of a missile oscillating in pitch at constant 
altitude with linearly varying aerodynamic coefficients: 

Iy 9 = C^q^Za + C m J^q, ( jL IQ + CmJ^q^jad (Ala) 


■where the not ah ions ( ) and ( ) refer ho the first and second derivatives 
with respect to time, respectively, and 


• - 
^ "51 


c mq 


= ^C m 



(Alb) 


C 







These 

range 


derivatives are assumed invariant with Mach number over the flight 
considered. From the diagram below. 
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e = a + 7 

and therefore 

e = d + 7 

•• .. •• 

e = a + 7 


(A 2 ) 


If Ct and C T . are neglected, the equation for rate of change of the 
Jj q. - 4 x 

flight direction is: 

mV 7 = 
or 

C I u Poo VAa 

7 = mV 2m 


(A 3 a) 

(A3b) 


Change variables from time to distance traveled, x, where the notations 
(»), (") refer to the first and second derivatives with respect to x, 

respectively. 


where 


so that 


7 = V 7 


2 « . t av 


7 = v 7 + 7 « 


dV _ _ c D p oo V A 

dt ~ m 2m 


• • <-y 1 1 

7 = V 2 7 - 


« C dPoo V 2 A7 


2m 


(A 4 a) 

(A 4 b) 

(A^c) 

(AAd) 


Substitute equation (A 3 b) into (A 4 a) 

. = 

1 2m 


(A 5 a) 


and 




(A 5 b) 


so that the derivatives of 7 with respect to distance are independent 
of velocity. Similarly, change variables from time to distance xn 
equation (A 2 ) . 
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^3 


e = a + 7 




0 = v 2 


tl 

a - 


(°D ' C 0 ^ “ " °^(^* 


a 


(A6) 


Substitute (A6) into (Al) and collect terms, 


tt 

a - 


_(% + Cl %) ' (% " CD ) 


P A' 

H oo 

2m 


[■ 


p Poo AZ . C Ia% f Pc/^ 2 

Cma l^T + V“ 


C D C Itt(^ 


Equation (A7) has the form: 

a + kja + = 0 

where k x and k 2 are constants defined by: 


a = 0 

(A7) 

(A8) 




C D - Ci^ + (^C mq + Cn y^jy- 


(A9) 


and 




The solution to equation (A 7 ) may be written in the following form: 


(All) 


where k x /2 is the damping constant, and C x and C 2 are determined by 
the amplitude and phase of the oscillation at x = 0 . The wave length 
and frequency of oscillation are given by 


a = e 


C x cos(k 2 - 



1/2 


x + C 2 sin(k 2 


ki 

h 


1/2 


X 


kh 
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2rt = 2rtf = 

A V 



(A12) 


Thus, A is constant if k x and k 2 are constant, but frequency of 
oscillation, f, is directly proportional to velocity. 

The effect of deceleration on the damping is given by the term, 
C^A/an in equation (A9) . The influence of deceleration on pitching 
wave length is shown by equation (A12) to depend on both ki and k 2 and 
thus on the additional terms, Cpp^/ 2m and Cp ) C]^(p cd A./2m) 2 . 


pitching motions are obtained by use of equations (A9) , (AlO) , and (A12) . 
They are: 


The equations for computing 





(A14) 
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APPENDIX B 

RELATIVE CONTRIBUTIONS OF AERODYNAMIC DAMPING AND 
DAMPING DUE TO DYNAMIC -PRESSURE VARIATION 


The significance of the pitch-damping coefficients measured and the 
variation of dynamic pressure due to density change on the flight charac- 
teristics of a missile as it descends from high altitude were investigated 
as described below. A full-scale missile (similar to the blunt 60° cone 
tested) was assumed, with the following geometric and aerodynamic 
properties: 

weight 125 slugs 

diameter ^ ft 

length 2.78 ft 


x cg 

l 


70 percent (from nose) 





constant, at values measured in the present investigations, 
dimen s i onle s s 


0 E entrance angle, 30° 

V E 20,000 ft/sec 

With the preceding conditions given, the following expression 
(ref. 20) may be used to describe the amplitude ratio of the missile in 
descending, power-off flight in the absence of aerodynamic damping: 

= (\ Wy* (Bl) 

cci VW 


The relations of reference 1 were used to describe the density varia' 
tion in the atmosphere and the velocity history of the missile; namely, 

c P p o A e -3y 
23m sin 


V = V E e 


(B2) 
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and 


where 


-By 

’o e = p > 

slugs/ ft 

, D = 0.0034 , 

slugs/ft 

B - 1 

ft" 1 

22,000 * 


Combining equations (Bl) and (B2) gives 


a max _ 
ad 


yi 


£ e -PYi 


- £ y £ e-py 

4 y e 2 


(B3a) 


where 

C D P 0 A 

C = - — 

2pm sin Q-g 


(B3t) 


Applying equation (B3) to the particular case described above leads to 
the curve labeled "dynamic-pressure effects, aerodynamic damping = 0" 
in figure 8 and represents the damping due to dynamic-pressure variation. 

For the purpose of comparing the rate of decay of pitching due to 
aerodynamic damping with that due to dynamic -pres sure change, it is use- 
ful to first consider level flight at constant speed at all altitudes 
within the range of interest. For flights of this type, the envelope 
of the pitching oscillations, a m ax> will he given by: 


<x max “ cx^e 2 (b4) 

where k x must be computed for Cp = 0 to avoid the introduction of 
dynamic-pressure-variation effects. The variable is a function of 

density alone since the aerodynamic coefficients have been assumed con- 
stant. The rate of change of amplitude with distance at every altitude 
is then 


1 ^%ax _ ki 

^max 2 


(B5) 
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where the transposition of o^ax to the denominator of the left-hand 
side gives the fractional change in a ma x per foot of travel. If it is 
assumed that the rate of decay of amplitude due to damping is at every 
altitude the same for descending flight as for level flight, then 
equation (B5) can he used to estimate the contribution of the aerodynamic 
damping to the pitching motion. The independent variable x must, how- 
ever, be replaced by the coordinate along the flight path y/sin 0g, 
where y is zero at the surface of the earth. 


sin dct-max _ kj 
a max dy 2 


From the definition of k x /2. 


ill 

2 



where 



-Dp 0 e _Py = -Dp 



(B 6 ) 


(B7a) 


(B7b) 


Combining equations (b 6) and (B7) and integrating 


gives 



^hnax 

^max 


P PQ 

sin 



p(p - Pj) k ij ~ k i 

a-max _ „ 3 sin 0 R _ _2(3 sin 0 R 

a i 


(b 8 ) 


(B9) 


Equation (B9) is plotted in figure 8 for the example missile. 

If equations (B3a) and (B9) are multiplied, the effect of both aero- 
dynamic damping and dynamic-pressure-variation damping should be obtained. 
The resulting equation: 
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total 


e 



(BIO) 


is also plotted in figure 8. Aerodynamic damping begins to affect the 
amplitude ratio at altitudes near 120,000 feet, and. becomes increasingly 
important at lower altitudes, so much so as to eliminate the diverging 
effect of dynamic-pressure variation at altitudes below 45,000 feet. 

While the present paper was in preparation, Mr. H. J. Allen solved 
equation (Al) for descending flight through the atmosphere, using the 
density and velocity variations given by equations (B2) , assuming 
constant aerodynamic coefficients, and neglecting only some second-order 
terms in the differential equation. The solution is given below, from 
reference (l6) , 



(Bll) 


where 



2 


(B12) 


and 



(B13) 


For large values of the Bessel function argument in equation (Bll) 
Allen found that the following equation was valid: 



a 




n/ it J k 3 + k 4 e 


4 


and, where the envelope value may be written: 


NACA RM A 56 IO 5 


k9 


^nax 


a E 


s ^3 e 


-By 


I — , 

4 It 4^-3 + k 4 e 4 


(B15) 


It was found that at the first half-cycle of oscillation, descending 
from infinite altitude, equation (B15) approximates equation (Bll) within 
1 percent, for the hypothetical missile chosen. The Bessel function argu- 
ment for this condition is approximately 3*83* The altitude corresponding 
to this condition is about 250,000 feet. 

In order to compare the results of Allen's solution with those pre- 
sented herein, <xi was chosen at y^ = 200,000 feet. Equation (B15) may 
then he written: 


a max _ je 

at 


-k 3 e 


-Byi 


-k 3 e - Py 

e J e 



(bi6) 


If aerodynamic damping is neglected in Allen's solution (^i.e., if 
and are equal to zero, or (c^ + c m^)(^) " C I U = 


equation (Bl6) reduces to equation (B3), Freidrich and Dore's expression. 
It can also he shown that if substitutions are made for c and D in 
in equation (BIO), this equation will he identical to Allen's solution 
in the form given in equation (Bl6) . 
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APPENDIX C 

METHOD OF CALCULATING BOUNDARY -LAYER THICKNESS 


The nomenclature used in this appendix is that given by Cohen and 
Reshotko in NACA TN's 3325 and 3326 (refs. 27 and 28). 

Equation (35) of reference 27 is, in transformed dimensionless 
variables, the general equation for velocity and temperature distribu- 
tion through the laminar boundary layer as a function of distance from 
the surface. When equation (35) is converted to the physical plane 
using dimensional variables, it becomes: 



where 

y normal coordinate 

6 momentum thickness 

t Q free- stream stagnation temperature 

t e local temperature outside boundary layer 

P pressure gradient parameter 

n correlation number 

S enthalpy function 

M e local Mach number outside boundary layer 

f» velocity ratio, JL ^°. 9 ?-?-..yeloclty in boundary layer 

^e local velocity outside boundary layer 

T| similarity variable 

As a particular value of y, the case of y = B can be considered. 

For this purpose, the boundary -layer thickness, 5 , is arbitrarily defined 
to agree with the definition used by Van Driest in his flat-plate calcu- 
lations (ref. 2 6); that is, y = 6 where u = 0.995 u e . Equation (Cl) 
may then be written as: 
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For a given value of 3 and Sy, the velocity ratio, f » , and enthalpy 
function, S, are tabulated as a function of the similarity variable, T), 
in reference 27* (The enthalpy function Sy is defined as: 

Sy = ("ty/to) - lj where by is the wall temperature.) The integral of 
equation (C2) is evaluated graphically by use of the tabulated data, 
thus leading to values of S /©. The procedure for determining the 
momentum thickness, 0, is clearly explained in reference 28. 

Experimental and theoretical pressure and velocity distributions 
were used to evaluate boundary- layer thickness, 5. For the hemisphere 
at a Mach number of 3*8, experimental values were used from reference 3. 
For the hemisphere at a Mach number of 8, simple impact theory was used 

(Cp = 2 sin 2 0 g ) . The shock-expansion method was used to obtain pressure 
and velocity distributions over the pointed ogive at a Mach number of 4. 

In order that the calculations for boundary- layer thickness be made 
with the minimum interpolation from the tables of reference 27. values 
of Sy used were those that were tabulated closest to the experimental 
conditions. Therefore, the calculations were made for models with 
slightly different heat-transfer conditions than were tested experi- 
mentally. The values used in the calculations and those that corre- 
sponded to experiment are shown in the table below, where the temperature 
ratios are given in terms of the notation of the present report; that is. 
Too is free-stream static temperature. 


Mach 

number 

Model 

Ty/Too 

(calc) 

T w/Too 

(exp) 

(calc) 

(exp) 

3.8 

hemisphere 

O .78 

1.0 

- 0.80 

-0.743 

4.0 

pointed ogive 

.84 

1.0 

-.80 

-.762 

8.0 

hemisphere 

2.76 

1.8 

-.80 

-.87 


Subsequent experience with the method of references 27 and 28 indi- 
cated that the effect on & of using Sy = -0.80 instead of Sy 
corresponding exactly to experiment would have been small. 
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(a) Round-nosed 60° cones. 


(h) Pointed 60° cone, hemisphere, and sphere. 


A-20728 



Figure 1.- Models. 
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1.600 

(Dia.) 


(c) Ogive, sine -wave -rounded ogive , and round -nosed 30° cone. 


Figure I. - Continued . 
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(d) Sine-wave-rounded ogive with Type III finish. 


Figure 1.- Concluded. 
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(a) Photomicrograph, Type I finish, 5 OOX. 



(b) Photomicrograph, Type III finish, ^>00X. 
Figure 2.- Photographs of model surfaces. 
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(c) Interferogram, Type I surface, 200X. 



(d) Interferogram, Type III surface, 200X. 


Figure 2.- Continued. 
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(e) Interferogram of stagnation pointy Type III surface, 200X. 

Figure 2.- Concluded. 


Angle of attack, a, degrees 




Distance, x, feet 


Figure 3.- Typical data for pitching oscillation of 


a round -nosed 60° cone. 
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O Round -nosed 60° cone 
□ Pointed 60° cone 
O Hemisphere 
A Sphere 
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Figure 4.- Drag data. 
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(a) Round -nosed cone, hemisphere, and sphere. 

Figure 5. - Comparison of computed and experimental total 

drag coefficients. 


Total drag coefficient, 
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(b) Pointed 60° cone and 60° cone - cylinder. 


Figure 5. - Concluded. 


Normal -force - curve slope, 
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O Experiment, round - nosed 60° cone 


Modified Newtonian theory 


Theory , 60° cone , reference 19 
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Figure 6. - Normal -force data. 



Center of pressure, 
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Experiment, round -nosed 60° cone 
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(Theoretical C M from reference 19) 
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Newtonian, round- nosed 60° cone, reference 18 
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Figure 7. - Center- of- pressure data. 
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Figure 8. - Relative contributions of dynamic - pressure vari- 
ation and aerodynamic damping during atmospheric 
entry of an example missile. 
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(a) Fully laminar flow off base, Moo = 



1 


(b) One side turbulent, = k. 

Figure 9 *- Shadowgraphs showing difference in flow field behind base with 

laminar and turbulent boundary layers. 
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(d) Smooth 60° cone- cylinder, = 8.1. 

Figure 9.- Concluded. 



* ^ 
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(a) Side-view, low sensitivity. 



(b) Top-view, high sensitivity. 

Figure 10.- Round-nosed 60° cone, Mo, = 8 . 3 . 
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(a) Dimensionless length of laminar run. 


O Visible turbulence 
X Burst waves 















( 

O 

X 

v 

9 ? 9 

o 

X 

o 

X 

o 

X 








X 











O 1 

-4 -3 -2 -I O I 2 3 4 


Sheltered Windward 

a R , degrees 

(b) Local Reynolds numbers at transition. 


Figure II.- Transition data from the round -nosed 60° cones with 
I surface finish at IVU = 8.3 and R^, « 4.8 x I0 6 . 
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(a) Dimensionless length of laminar run. 
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(b) Local Reynolds numbers at transition. 


Figure 13.- Transition data, round-nosed 60° cones at Mod =4.0 


and Ra>= 3.8 x ICf. 
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(a) Smooth. 



(b) Roughened. 


Figure 14 .- Pointed 60° cone. Moo = 8.3. 
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(b) Roughened, Moo = 8.2. 


Figure 15 


Hemispheres 
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Figure 15.- Concluded. 
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(a) Normal shadowgraph. 



(b) Less-sensitive shadowgraph. 
Figure 1 6 .- Ogive, = b. 
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Figure 17. - Transition data from pointed ogives. 
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Figure 18 .- Sine -wave -rounded ogive 
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Figure 19.- Length of laminar runs on sine-wave-rounded 
ogives compared with those on pointed ogives. 
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Figure 20. - Comparison of local transition Reynolds numbers on sine-wave- rounded 

and pointed ogives; 4, 4 x 10? oo 
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(a) Small angle of attack. 



Figure 21.- Round-nosed 30° cone. Moo - 


1 . 0 - 


<5 

t- 

w 


0^- 

-20 


1.5 




I 


? 


-10 


0 io 

a , degrees 
(a) Right 


Sheltered 


Windward 


a , degrees 

(c) Top 


20 







( 

0 

X 

6 

9 

. X 


1 

2 

X 

X 

k 





O Visible turbulence 
X Burst waves 
| Eddies in shear 

0 


^ layer 




0 


"1 

( 

J 

: 



-20 


Angle of attack 


-10 


0 10 
a , degrees 
(b) Left 


Sheltered Windward 

a , degrees 
(d) Bottom 


20 


| 





0 

0 c 

5 1 

/-s 



XX 

u 

X 



03 

VJl 


Figure 22.- Length of laminar run on round-nosed 30° cone 
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Figure 23. - Variation of roughness parameter over surface 

of pointed ogive. 
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Figure 24.- Variation of critical -roughness parameter with 
boundary -layer -edge Mach number for pointed bodies. 



Position on surface, s /d 


Figure 25 .- Boundary -layer thickness development on spherical nose. 
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Figure 26- Pressure - distribution data for the round -nosed 80° cone of reference 9. 
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(a) Leading-edge thickness = 0.0003 inch. 



(h) Leading-edge thickness = 0.0100 inch. 


Figure 27.- Comparison of turbulent boundary-layer initial growth rates with and without an appre- 
ciable inviscid shear layer due to leading-edge bluntness; Moo = Roo = 17 X 10 6 , based on 

cylinder length. 
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Figure 28. -Effect 


of (CoA/m) on Reynolds numbers 
entering the atmosphere. 
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